Driving particle current through narrow channels using classical pump 
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We study a symmetric exclusion process in which the hopping rates at two chosen adjacent sites 
vary periodically in time and have a relative phase difference. This mimics a colloidal suspension 
subjected to external time dependent modulation of the local chemical potential. The two special 
sites act as a classical pump by generating an oscillatory current with a nonzero T>C value whose 
direction depends on the applied phase difference. We analyze various features in this model through 
simulations and obtain an expression for the T>C current via a novel perturbative treatment. 
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A number of studies have used the idea of Thouless 
adiabatic pumping T| to generate T>C current responses 
to AC driving fields in quantum systems. This mech- 
anism has been used experimentally to generate charge 
and spin 0] current in systems such as quantum dots 
and nanotubes and studied theoretically in various quan- 
tum systems without Q and with [f| interactions. One 
may ask if this principle can be used to drive current in 
classical systems. It was shown in [(| that inelastic scat- 
tering can enhance pumping which suggests that similar 
effects can be obtained in purely classical systems. In- 
deed classical pumping of particles and heat has been 
achieved in mesoscopic systems with large inelastic scat- 
tering [61, stochastic models |7y8|,|9( and seen in experi- 
ments [10| . As discussed in II ll these pump models are 
similar to Brownian ratchets [12j | where non-interacting 
particles placed in spatially asymmetric potentials that 
vary periodically in time and acted upon by noise exe- 
cute directed motion. However for the model described 
below, particle interactions seem necessary for pumping. 

Here we show how the designing principle of quantum 
pumps may be applied to drive classical particles such 
as micron-sized charged colloidal particles confined in a 
closed narrow tube. Systems of colloidal particles with 
excluded volume interactions diffusing in a narrow chan- 
nel have been modeled in the recent past [l3|, [3] by one- 
dimensional symmetric exclusion process (SEP) in which 
hard-core particles attempt to hop to an empty neigh- 
bor with equal rate [l5j |. Recently much attention has 
been given to nonequilibrium steady states of driven SEP 
lil 17 1, in which particles can enter or leave the system 
at the boundaries, to study large deviation functional [18| 
and current fluctuations (l9( in nonequilibrium systems. 

In this Letter, we introduce a SEP in which the hop- 
out rates at some lattice sites are chosen to be time- 
dependent with a relative phase difference between dif- 
ferent points. At a coarse-grained level, this choice mod- 
els oscillating voltages applied to the colloidal system. 
We find that an oscillatory particle current is generated 
across a bond which interestingly has a nonzero T>C value. 
Unlike driven SEP, our model conserves the total number 
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FIG. 1: (Color online) T>C density profile pi across the ring 
for f = 0.3, fx = 0.2, u = 0.2n and cj> = ir/2 at half filling 
for two system sizes. Inset: T>C current J ~ l/L as shown by 
solid line of slope — 1. 



of particles and due to the periodically varying hopping 
rates, the system reaches in the long time limit a periodic 
time-dependent state. 

Here we will mainly focus on the model with a single 
pump which we define as consisting of two adjacent sites 
where the hopping rates vary periodically and with a 
phase difference. The hopping rate at all other sites is 
constant. Our Monte-Carlo simulations indicate that the 
magnitude of the T>C current J depends sinusoidally on 
the phase difference between the rates at these two sites, 
and varies non-monotonically with the frequency of the 
drive vanishing in both zero and infinite frequency limits. 
We also obtain an analytical understanding of the system 
by developing a perturbation theory in the strength of 
the time-dependent part of the hopping rates, similar to 
(20| . We show that the leading order contribution to J 
is obtained at second order in this perturbing parameter, 
and find an explicit expression for it (see Eq. (|15p ) which 
captures various features seen in the simulations. 

The model is defined on a ring with L sites. A site 
I = 1, 2..., L can be occupied by ni — or 1 particle, and 
the system contains a total of N = pL particles where p 



it/2 n 3x12 2it 



FIG. 2: (Color online) Plot of J vs. with other parame- 
ters as in Fig. [T] The solid lines are from the perturbation 
theory.Inset shows plot of J vs. to 



is the density. A particle at site / hops to an empty site 
either on the left or right with rate ui = fo for alU ^ 1, L 
and ui = f + fx sin(wt), u L = / + /i sin(wt + </>). For 
fx = 0, the above model reduces to the SEP with peri- 
odic boundary conditions, many of whose properties are 
known exactly [lij]. The steady state obeys the equilib- 
rium condition of detailed balance so that the current 
across a bond is zero. Due to this, the steady state 
measure is uniform and the fc-point correlation function 

C£l,...,i k = (n h n h ...n lh ) = (£:*)/(£)• As we will de- 
scribe later, this knowledge of the exact steady state in 
the absence of time-dependent term in ui allows us to set 
up a perturbation expansion in fx of various observables. 

We first discuss the results of our Monte-Carlo simula- 
tions of the model defined above for generic values of the 
driving frequency u> and the phase <$>. We first look at the 
VC values of the current and local densities. In Fig.Q] we 
have plotted the local particle density profile for system 
sizes L = 16 and 64. We find a linear density profile in 
the bulk of the system with a discontinuity at the pump- 
ing sites. The linear density profile implies a VC current 
in the system. The inset in the figure shows that this 
current scales as 1 /L as in driven SEP IH, . In Fig. [2] 
we plot the current as a function of the applied phase dif- 
ference for different system sizes. The prediction of the 
perturbation theory (described below) is also shown for 
comparison. As in the case of quantum pumps, we find a 
sinusoidal dependence of the current on the phase </>. The 
inset of Fig. [2] shows the dependence of the current on the 
driving frequency to again for L = 16, 64. At lo = 0, the 
system is in equilibrium and there is no current [2l( . In 
the opposite limit lo — > oo, the system does not have time 
to react to the rapidly changing fields and we may again 
expect zero current. For any finite lo, the VC current is 
nonzero and behaves as 




FIG. 3: (Color online) Plot of time-dependent densities at 
the four sites of a L = 4 lattice. In the initial configuration, 
sites 1 and 2 have one particle each and other sites are empty. 
The averages over one time period give: px = 0.503493, p2 = 
0.498702, p 3 = 0.497417, p 4 = 0.500388 and J = 0.000514. 
The points show the curve p + fxPx + fxPx • [Parameters: 
fo = 0.4, fx = 0.1, 4> = tt/2 and lo = 0.2tt]. 



where lo* is the frequency at which the current peaks. A 
similar non-monotonic behavior in current has been seen 
in experiments on ion pumps HI . Since J ~ lo as lo — > 0, 
this means that a finite number of particles are circulated 
even in the adiabatic limit [l(| • This suggests a perturba- 
tion expansion in lo and will be discussed elsewhere. We 
also find that lo* scales linearly with / i is independent 
of L and hence of the typical relaxation time r = L 2 / f 
of the system. 

Till now we have looked at various time averaged quan- 
tities. It is also of interest to look at the time depen- 
dences. In Fig. [3] we show the ensemble averaged local 
densities as a function of time obtained by an exact nu- 
merical integration of the L = 4, N = 2 master equation. 
We find that the local densities oscillate about their re- 
spective VC value ( different from p = 0.5) at large times, 
and a nonzero VC current circulates around the ring. Al- 
though the local VC density is close to p, the amplitude 
of the AC part is quite large. We also show a comparison 
of the exact density pi (t) with the approximate density 
P + fiPt obtained from the perturbation theory 

(see Eqs. © and (HI])). We find that the AC density con- 
tains the driving frequency and also higher harmonics. 

We now proceed to describe a perturbation theory 
within which we will calculate an expression for the VC 
current J in the bulk of the system. Let us denote the 
site hopping rates by 



Uj 



fo + fm =f + fi(u ie iut + vle- wt ) , (2) 



LO < LO*, 
LO > LO* 



(1) 



where the site-dependent complex amplitude v\ — 
—ill, Vl = —ie l ^/2 and all other Ui's are zero. We con- 
sider perturbative expansions of various quantities of in- 
terest with f\ as the perturbation parameter about the 
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homogeneous steady state corresponding to /i = 0. Thus 
we have 



Pl {t) = (n,(t)) = P + X;/M fc) (*) 

k=l 

OO 

CU(t) = (m(t)M*)) = cS + E/i fc c£(*) J (3) 

fe=i 

and similar expressions for higher correlations. The cur- 
rent in a bond connecting sites I and Z + 1 is given by 
Ji,i+i(t) = ui(pi-Ci t i + i)-ui + i(pi+i-Ci t i + i) which can 
also be expanded in a perturbation series using Eq. ([3]). 
Averaging Jij + i(t) over a time period for I ^ L — 1, 1, L, 
we obtain constant VC current in the bulk [221, 



fc=i 
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where pf" is the time-averaged density at the k th order. 
To obtain J to leading order in /i, we start with the 
exact equation of motion for the density pi (23l | 

<9pz 

— + 2Uip/ - - U i+ ip ; + i 

OT 

= ui(Ci-x,i + Cij+i) - ui+iQj+i - ui-iCi-ij . (5) 

Then the density p\ ' at the fc th order evolves according 
to 



9Pi 



(fc) 



dt 



(fc-i) 



(fc-i) 



where A;(^ = — 2<?; + <?;-i defines the discrete Lapla- 
cian operator. Thus, the density at the k th order is ob- 
tainable in terms of the density and the two-point corre- 
lation function at (k — l) th order. 

At first order, the above equation simplifies to give 



dt 



vanish (using the condition J2i Pi (0) — 0) while the 
particular integral has the asymptotic form 



pf } (t) 



with the amplitudes given by 



i . i 
Q-+z + + a_z_ 



I = 2,3...,L- 1 



(8) 
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where z + = y/2 + [(y/2) 2 - l] 1 / 2 , z_ = l/z + , y = 2 + 
iw// and a± = (~ir uj/f§ ){v x z T + v L )/[{z± - 1)(% ~ 
z±)]. The amplitude of the density oscillations decay 
exponentially with distance from the pumping sites. 

We also need the two-point correlation function C|„ 
at C(/i) appearing in Eq. ^ for k = 2. Due to hard- 
core interaction, the correlation function C; !?n for nearest 
neighbors I and m obeys a different equation of motion 
when \l — m\ ^ 1. Taking care of this complication, it 
is straightforward to write down the corresponding evo- 
lution equations [23} and carry out the perturbative ex- 
pansion. At C(/i) in perturbation theory, we then get 
(for 1 < I < L, I + 1 < m < L) 



dC, 



(i) 



dt 



/o(A, + Ajcf] = k (A lVl + A, 



dC, 
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ot 
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1,1+1 



^1-1,1+1 



_r»(i) 
^1,1+2 

Vl+2 - Vl 



vi+i) (9) 



where fco = C 



(0) 
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The computation of even 



the homogenous solution of the above set of equations 
is in general a nontrivial task because of the form of the 
equations involving nearest neighbor indices and requires 
a Bethe ansatz or dynamic product ansatz How- 
ever the long time solution can still be found exactly and 
is given by: 



f A lP [ 1] = r A 



(7) 



a 



/.m(*) = ^-iPl 



(1) 



p -iu>t 



where r = p — C t 
both sides of Eq. 



(0) 



Pl+L an< ^ density conservation, we find p[^ 



(10) 

In the steady state, integrating where A^ n — (fc /r )(A| 1 ^ + Am)- It is easily verified 

that this satisfies Eq. (0 for all I, m. Interestingly, if one 
were to assume that the two-point function decouples, 
then a perturbation expansion to order /1 will also yield 
the preceding equation in the thermodynamic limit. To 
determine whether the system indeed has a product mea- 
sure requires a more detailed analysis of the higher order 
terms in the perturbation series and higher correlations, 
and will be discussed elsewhere. 

We are now in a position to find the density profile 
p[ 2) (t). Using Eqs. © and (HUJ) in Eq. © for k = 2, we 
get 



over a time period, and using p\ 

0. Hence 

at C(/i), no VC current is generated. This is consistent 
with our expectation from linear response theory, namely 
that an AC field gives an AC response. Thus, we need to 
consider the next order term p\ (t) in the perturbation 
expansion of the density to obtain current. As Eq. © 
for k = 2 involves density at first order, we first solve for 
Pl (t) in Eq. 0. The solution of this inhomogeneous 
linear equation is a sum of a homogeneous part which 
depends on initial conditions, and a particular integral. 
At long times, the homogeneous part can be shown to 



pf\t)=p?HAf>^ 



A^e 



■i2ut 



(11) 
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where the DC part p\ ' obeys the equation 

-ut(4Xi + 4V+i) + ^-iA a) i, + ^ + i4V + i](i2) 

Noting that the right hand side of the above equation is 
zero for most I, we get the following solution 

pf ] = sl + h, Z = 2,...,L-1 

Pi 2) = s + ft + ^ReKCAg-^)] 
Jo 

= s l + h|- Rc K(4\l-4 1, )]> (13) 

Jo 

where the slope s of the linear density profile is given by 

s = j^r M^l-A^ L ) + vt{A%-A^_ ltL )] , (14) 
Wo 

and the intercept h can be found using the particle con- 
servation condition p\ =0. 

Our final result for the leading order contribution to 
the VC current in Eq. (jJJ is J = —fof%s and, using the 
form of A^L, can be written as 



./ = - f £ 

Jo 



koto sin cj) 
2L 



,L-1 



Rc 



1 



out) 



(15) 

We observe that the 1/L dependence of the VC current 
and its sinusoidal variation with the phase are captured 
by the perturbation theory at 0{f1). One can also check 
from the above expression that the frequency dependence 
of current in Eq. ([1} is obtained with the typical fre- 
quency oj* w 2/o, I>1. As shown in Fig. ([2]) the above 
perturbative solution is already in a very good quantita- 
tive agreement with the simulation results. 

To summarise, we have studied a stochastic model of 
hard-core particles on a ring in which periodically vary- 
ing hopping rates can induce a VC current. A variation of 
our model defined on an open chain connected to particle 
reservoirs also shows nonzero VC current. Simulations 
and a novel perturbative analysis presented here show 
that many of the qualitative features are similar to that 
seen in quantum pumps. However unlike in quantum 
pumps, interactions play an essential role in our model. 
For non-interacting particles, noting that the bond cur- 
rent Ju+i = uipi — u/_|_ipz_|_i summed over all bonds is 
zero and using the uniformity of J, it follows that no VC 
current is generated in this case. Our model can be gen- 
eralised to include several pumps and pumps consisting of 
several sites, and these features could increase pumping 
efficiency. 
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A number of studies have used the idea of Thouless 
adiabatic pumping lj to generate DC current responses 
to AC driving fields in quantum systems. This mecha- 
nism has been used experimentally to generate charge 
and spin [3| current in systems such as quantum dots and 
nanotubes and studied theoretically in various quantum 
systems without Q and with [B[ interactions. One may 
ask if this principle can be used to drive current in classi- 
cal systems. It was shown in |6| that inelastic scattering 
can enhance pumping which suggests that similar effects 
can be obtained in purely classical systems. Indeed clas- 
sical pumping of particles and heat has been achieved 
in certain stochastic models 0JH, Q and seen in experi- 
ments [Ioj |, As discussed in these pump models are 
similar to Brownian ratchets [12| where non-interacting 
particles placed in spatially asymmetric potentials that 
vary periodically in time and acted upon by noise exe- 
cute directed motion. However for the model described 
below, particle interactions seem necessary for pumping. 

Here we show how the designing principle of quantum 
pumps may be applied to drive classical particles such 
as micron-sized charged colloidal particles confined in a 
closed narrow tube. Systems of colloidal particles with 
excluded volume interactions diffusing in a narrow chan- 
nel have been modeled in the recent past [H, Q3] by one- 
dimensional symmetric exclusion process (SEP) in which 
hard-core particles attempt to hop to an empty neigh- 
bor with equal rate [HI]. Recently much attention has 
been given to nonequilibrium steady states of driven SEP 
HI 17 1, in which particles can enter or leave the system 
at the boundaries, to study large deviation functional [l8| 
and current fluctuations [191 ] in nonequilibrium systems. 

In this Letter, we introduce a SEP in which the hop- 
out rates at some lattice sites are chosen to be time- 
dependent with a relative phase difference between dif- 
ferent points. At a coarse-grained level, this choice mod- 
els oscillating voltages applied to the colloidal system. 
We find that an oscillatory particle current is generated 
across a bond which interestingly has a nonzero DC value. 
Unlike driven SEP, our model conserves the total number 
of particles and due to the periodically varying hopping 
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FIG. 1: (Color online) T>C density profile pi across the ring 
for f = 0.3, fx = 0.2, u = 0.2n and cj> = ir/2 at half filling 
for two system sizes. Inset: T>C current J ~ l/L as shown by 
solid line of slope — 1. 



rates, the system reaches in the long time limit a periodic 
time-dependent state. 

Here we will mainly focus on the model with a single 
pump which we define as consisting of two adjacent sites 
where the hopping rates vary periodically and with a 
phase difference. The hopping rate at all other sites is 
constant. Our Monte-Carlo simulations indicate that the 
magnitude of the DC current J depends sinusoidally on 
the phase difference between the rates at these two sites, 
and varies non-monotonically with the frequency of the 
drive vanishing in both zero and infinite frequency limits. 
We also obtain an analytical understanding of the system 
by developing a perturbation theory in the strength of 
the time-dependent part of the hopping rates, similar to 
(20| . We show that the leading order contribution to J 
is obtained at second order in this perturbing parameter, 
and find an explicit expression for it (see Eq. (fl"5|) ) which 
captures various features seen in the simulations. 

The model is defined on a ring with L sites. A site 
I = 1, 2..., L can be occupied by ni = or 1 particle, and 
the system contains a total of N = pL particles where p 
is the density. A particle at site / hops to an empty site 



it/2 n 3x12 2it 



FIG. 2: (Color online) Plot of J vs. <j!> (main) and J vs. w 
(inset) with other parameters as in Fig. [T] The solid lines are 
from the perturbation theory (see Eq. I|15p ). 



either on the left or right with rate ui = fo for alU ^ 1, L 
and mi = f + fi sm(u>t), u L = f + fx sin(wt + cf>). For 
fx = 0, the above model reduces to the SEP with peri- 
odic boundary conditions, many of whose properties are 
known exactly [l5(. The steady state obeys the equilib- 
rium condition of detailed balance so that the current 
across a bond is zero. Due to this, the steady state 
measure is uniform and the fc-point correlation function 

C h} 2 ,-,l k = (n h n h ...n lh ) = As we wil1 de- 

scribe later, this knowledge of the exact steady state in 
the absence of time-dependent term in ui allows us to set 
up a perturbation expansion in fx of various observables. 

We first discuss the results of our Monte-Carlo simula- 
tions of the model defined above for generic values of the 
driving frequency u> and the phase 0. In Fig. [TJ we have 
plotted the local particle density profile for system sizes 
L = 16 and 64. We find a linear density profile in the 
bulk of the system with a discontinuity at the pumping 
sites. The linear density profile implies a T>C current in 
the system. The inset in the figure shows that this cur- 
rent scales as 1/L as in driven SEP [H, [ll] . In Fig. [2] we 
plot the current as a function of the applied phase dif- 
ference for different system sizes. The prediction of the 
perturbation theory (described below) is also shown for 
comparison. As in the case of quantum pumps, we find a 
sinusoidal dependence of the current on the phase <f>. The 
inset of Fig. [2] shows the dependence of the current on the 
driving frequency lo again for L = 16, 64. At lo = 0, the 



system is in equilibrium and there is no current 2l|. In 
the opposite limit lo — > oo, the system does not have time 
to react to the rapidly changing fields and we may again 
expect zero current. For any finite lo, the T>C current is 
nonzero and behaves as 



J 



LO < OT 



(1) 



1/lO , LO > LO* , 

where lo* is the frequency at which the current peaks. A 



FIG. 3: (Color online) Plot of time-dependent densities at 
the four sites of a L = 4 lattice. In the initial configuration, 
sites 1 and 2 have one particle each and other sites are empty. 
The averages over one time period give: pi = 0.503493, p2 = 
0.498702, p 3 = 0.497417, p 4 = 0.500388 and J = 0.000514. 
The points show the curve p + fxPx + fxPx ■ [Parameters: 
fo = 0.4, fx = 0.1, 4> = tt/2 and lo = 0.2tt]. 



similar non-monotonic behavior in current has been seen 
in experiments on ion pumps [10(. Since J ~ lo as to — > 0, 
this means that a finite number of particles are circulated 
even in the adiabatic limit [l(| • This suggests a perturba- 
tion expansion in lo and will be discussed elsewhere. We 
also find that lo* scales linearly with /o, is independent 
of L and hence of the typical relaxation time r = L 2 / fo 
of the system. 

Till now we have looked at various time averaged quan- 
tities. In Fig. [3] we show the ensemble averaged local 
densities as a function of time obtained by an exact nu- 
merical integration of the L = 4, N = 2 master equation. 
We find that the local densities oscillate about their re- 
spective T>C value ( different from p = 0.5) at large times, 
and a nonzero T>C current circulates around the ring. Al- 
though the local VC density is close to p, the amplitude 
of the AC part is quite large. We also show a comparison 
of the exact density pi (i) with the approximate density 
p + fiPl obtained from the perturbation theory 

(see Eqs. © and ([IT])). We find that the AC density con- 
tains the driving frequency and also higher harmonics. 

We now proceed to describe a perturbation theory 
within which we will calculate an expression for the T>C 
current J in the bulk of the system. Let us denote the 
site hopping rates by 



ui = fo + fiV l = f + f 1 (p l e iuJt + , 



') , (2) 



where the site-dependent complex amplitude v\ — 
—i/2, vl = — and all other vi's are zero. We con- 
sider perturbative expansions of various quantities of in- 
terest with f\ as the perturbation parameter about the 
homogeneous steady state corresponding to fx — 0. Thus 



3 



we have 

oo 

Pl (t) = (n l (t)) = p + J2fiPl k \t) 



fe=i 



(ni(t)n m (t)) = C[% + V /*C<*> (t) , (3) 



fe=i 



and similar expressions for higher correlations. The cur- 
rent in a bond connecting sites I and I + 1 is given by 
Ji,i+i(t) = ui(pi - Cu+i) - u l+1 (pi +1 - Ci,j + i) which can 
also be expanded in a perturbation series using Eq. ([3]). 
Averaging Ji t i + x{t) over a time period for £ ^ Z — 1, 1,£, 
we obtain constant T>C current in the bulk, 



(4) 



fe=i 



where p| fc ^ is the time-averaged density at the k th order. 
To obtain J to leading order in fx, we start with the 
exact equation of motion for the density pi [2^ | 



9* 



+ tuipi - ui-ipi-i - ui +1 pi + i 



= ui(Ci-ij + Cu+i) — ui+iCij+i - ui-iCi-ij .(5) 

(k) 

Then p\ evolves according to 

/oA/p, + 2v i p l - Vi^xPi-x - v i+iPi+i 

- + + 1 ) ^ v i-i c i-x,i -vi+i^i,i+x >W 

where Aigi = gi + x — 2gi+gi_x defines the discrete Lapla- 
cian operator. Thus, the density at the k th order is ob- 
tainable in terms of the density and the two-point corre- 
lation function at (k — l) th order. 

At first order, the above equation simplifies to give 



^-/oA^ (1) =r A^ , 



(7) 



where ro = p — C\ j . In the steady state, integrating 
both sides of Eq. (0 over a time period, and using p"p = 
Pl+L and density conservation, we find pj = 0. Hence 
at C(/i), no VC current is generated. This is consistent 
with our expectation from linear response theory, namely 
that an AC field gives an AC response. Thus, we need to 
consider the next order term p\ (t) in the perturbation 
expansion of the density to obtain current. As Eq. ([6|) 
for k — 2 involves density at first order, we first solve for 
Pl(t) in Eq. 0. The solution of this inhomogeneous 
linear equation is a sum of a homogeneous part which 
depends on initial conditions, and a particular integral. 
At long times, the homogeneous part can be shown to 
vanish (using the condition pj 1 ^ (0) = 0) while the 
particular integral has the asymptotic form 



pf\t) 



(8) 



with the amplitudes given by 



(i) _ 
I ~ 

(i) _ 



l 

A<P = 



cx+z + + a_z_ 



a + z+ + a_z^ — 



I = 2,3..., L- 1 
rpvx 
fo 

fo 



where z + = y/2 + [(y/2) 2 - l] 1 / 2 , z_ = l/z+, y = 2 + 
iu}/f a and a± = (-ir uj/ f$)(v 1 z T + u L )/[(z^ ~ l)(z T ~ 
z±)]. The amplitude of the density oscillations decay 
exponentially with distance from the pumping sites. 

We also need the two-point correlation function cPjL 
at C(/i) appearing in Eq. © for k — 2. Due to hard- 
core interaction, the correlation function C^ m for nearest 
neighbors I and m obeys a different equation of motion 
when |Z — m\ ^ 1. Taking care of this complication, it 
is straightforward to write down the corresponding evo- 
lution equations [22| and carry out the perturbative ex- 
pansion. At O(fx) in perturbation theory, we then get 
(for 1 < I < L, I + 1 < m < L) 



dC, 



(i) 



dt 



/o(Aj + A m )C ! (1 ] = k (A m + A m v m ) , 



dt 



- + Jo y^ij+x ^l-X,l+X ^1,1+2 ) 

= k (vi_x + vi+2 -vi- vi+i) (9) 



where kg 



c (o) _ c (o) 



The computation of even 



the homogenous solution of the above set of equations 
is in general a nontrivial task because of the form of the 
equations involving nearest neighbor indices and requires 
a Bethe ansatz or dynamic product ansatz [l7| • How- 
ever the long time solution can still be found exactly and 
is given by: 



(1) -iwi 



(10) 

where = (fcoAoX^; 1 ^ + A^m). It is easily verified 

that this satisfies Eq. © for all /, m. Interestingly, if one 
were to assume that the two-point function decouples, 
then a perturbation expansion to order fx will also yield 
the preceding equation in the thermodynamic limit. To 
determine whether the system indeed has a product mea- 
sure requires a more detailed analysis of the higher order 
terms in the perturbation series and higher correlations, 
and will be discussed elsewhere. 

("2) 

We are now in a position to find p\ (£). Using Eqs. ((HJ) 
and (TO) in Eq. © for k = 2, we get 



pf\t) = pf ) + A^e^ + 4 (2 V« (11) 

-(2) 

where the VC part p\ obeys the equation 



f = 2 Re[2^ - vUA?\ ~ 



foAtp^ = 2 Re{2v*A\ L> - v^AQ - v, .1 



4 



Noting that the right hand side of the above equation is 
zero for most I, we get the following solution 
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pf ] = sl + h , 1 = 2,...,L-1 



p[ 2) = s + h + jr Re[v^(A\^ - A\ L) )} 
Jo 

P ( l } = S L + h+^R e y L (A^ L ~A^)}, (13) 
Jo 

where the slope s of the linear density profile is given by 

s = M^(A^l-A^) + ^ L (A^ L -A^_ 1>L )] , (14) 
Wo 

and the intercept h can be found using the particle con- 
servation condition J2i Pi = 0- 

Our final result for the leading order contribution to 
the VC current in Eq. ^ is J — —foffs and, using the 
form of A^ v can be written as 



(i) aV)\ 



koto sin (j) 
2L 



Re 



1-zi 



(15) 

We observe that the 1/L dependence of the DC current 
and its sinusoidal variation with the phase are captured 
by the perturbation theory at O(ff). One can also check 
from the above expression that the frequency dependence 
of current in Eq. (flj is obtained with the typical fre- 
quency uj* « 2/ , L>1. As shown in Fig. © the above 
perturbative solution is already in a very good quantita- 
tive agreement with the simulation results. 

To summarise, we have studied a stochastic model of 
hard-core particles on a ring in which periodically vary- 
ing hopping rates can induce a T>C current. A variant of 
our model defined on an open chain connected to particle 
reservoirs also shows nonzero T>C current. Simulations 
and a novel perturbative analysis presented here show 
that many of the qualitative features are similar to that 
seen in quantum pumps. However unlike in quantum 
pumps, interactions play an essential role in our model. 
For non-interacting particles, noting that the bond cur- 
rent Ju+i = uipi — ui + ipi + i summed over all bonds is 
zero and using the uniformity of J, it follows that no DC 
current is generated in this case. Our model can be gen- 
eralised to include several pumps and pumps consisting of 
several sites, and these features could increase pumping 
efficiency. 

K.J. thanks A. Rakos for useful discussions, Raman 
Research Institute for kind hospitality and Israel Sci- 
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